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In the framework of the generalized mean field theory, conditions for arising the ferromagnetic 
state in a two-dimensional diluted magnetic semiconductor and the features of that state are 
defined. RKKY-interaction of magnetic impurities is supposed. The spatial disorder of their 
arrangement and temperature alteration of the carrier degeneracy are taken into account. 



Diluted magnetic semiconductors, such as Gai_ x Mn.2; As , are broadly investigated in con- 
nection with their potential for new electronics developments and, especially, spintronics. For 
explanation of the ferromagnetism in those compounds the known RKKY-mechanism of the 
indirect exchange interaction is recruited PjJ |2] which leads to the correct estimating the Curie 
temperature in the framework of the traditional mean field theory. Mn-atoms (with concentra- 
tion Nfj) substituting for Ga-atoms introduce in the system the own magnetic moments and, in 
addition, as acceptors deliver free holes (with concentration n). It is precisely those holes be- 
come to be carriers responsible for the interaction. However, the equality of the concentrations 
n = keeps only at low Mn-concentrations (x < 0.05) PQI2], so that the carrier concentration 
is usually less than the concentration of magnetic impurities: n = jN^ where the coefficient of 
the impurity efficiency 7 < 1 and falls 1 with increasing N^. 

Nevertheless, the concentration of magnetic impurities, delivering carriers, in actual systems 
is usually so high that the impurity band is formed which at x > 0.01 merges into the valence 
band 5J. Even though, the carrier concentration occurs to be not so high that one could consider 
them as highly degenerated ones within the whole (being of interest) range of relatively high 
temperatures. Furthermore, it is important that in that range the carrier concentration is 
almost independent of the temperature: n = 7^ m Const. 

Although two-dimensional structures represent the most natural systems for the embedding 
in the traditional semiconductor technology, almost all theoretical works are dealt with the 
three-dimensional systems of degenerate carriers. Similarly, most of experimentally studied 
systems are three-dimensional ones. 

The objective of the present paper is considering magnetic features of two-dimensional 
semiconductor systems with magnetic impurities interacting by RKKY-mechanism via carriers 
of arbitrary degeneracy. That problem has been recently considered in the paper [6J where it 
has been shown that reducing the system dimension (from 3D to 2D) results in the significant 
lowering of the Curie temperature (under equivalent parameters). We think though that there 
are some inexactitudes in the paper. Firstly, authors has neglected of the temperature 
dependence £f{T) of the carrier Fermi energy that is inadmissible under the intermediate 
degeneracy when e F jk B T ~ 1 (and, all the more, at Bp < 0). In addition, the disorder in 
the arrangement of magnetic impurities has been considered in the framework of the mean 



^^One could control the relative hole concentration (i.e., the 7- value) by simultaneous introducing non- 
magnetic acceptors (for instance, Be 0E1) or choosing the temperature of the film growth 0. 



field theory (continual or the so called lattice one) where the averaged impurity configuration 
is ordered one and the scattering of local effective magnetic fields is not considered. As we 
demonstrate below, taking into account the mentioned things influences the results significantly. 

The general expression for RKKY interaction in two-dimensional systems is known |B1 E| , 
however the relevant results refer to the degenerate carriers and have to be generalized for the 
case of the intermediate degeneracy \ep/kBT\ ~ 1 {ep is the carrier Fermi energy). 

In a few papers [HH ITT] (and later in [T2|) the systems with non-degenerate carriers have been 
considered and it has been shown that the energy of RKKY interaction in the non-degenerate 
case is exponentially falls with the distance (in contrast to the degenerate case where the fall is 
power behaved): J oc exp(— r 2 /A|), where A T = h/(2mk B T) 1 / 2 is the mean (thermal) de Broglie 
carrier wave length. However, those papers, first, refer to the three-dimensional systems and, 
second, they lack of taking into account the random arrangement of magnetic impurities. 

It is known that the traditional mean field theory does not provide the adequate description 
of a disordered (random) system of magnetic moments In the present paper, we shall 

use the generalized mean field theory ^3] for systems with the indirect interaction of magnetic 
impurities taking into account the randomness of their spatial arrangement. We shall use 
Ising approximation and suppose that the indirect coupling between magnetic moments of 
impurity atoms is realized by means of RKKY interaction which is replaced by the effective 
magnetic field, whereupon system properties are described with the help of the distribution 
function of local values of the field arising as a result of magnetic ions' coupling with their 
own surroundings. In real systems, the scattering of those fields proves to be so substantial 
that RKKY interaction makes the magnetic ordering possible at lower temperatures only (as 
compared to those predicted by the traditional mean field theory). Dependencies of magnetic 
properties on the carrier concentration turn out to be unusual, as well. 

In (ZllHlini, the expression has been derived for the energy w(r) of indirect RKKY interaction 
for two parallel spins Si, S2 of magnetic ions spaced at the distance r in the two-dimensional 
system with degenerated carriers: 

2 kp 
w{r) = (^\ F(r) S X S 2 , F{r) = - J kN (kr)J Q (kr)dk, (1) 



where J ex is the exchange energy for interaction of a spin with a free charge carrier of the mass 
m, N is the concentration of lattice atoms (N = 1/a 2 for the square lattice of the period a); J^i 
iVo,i are Bessel functions. To generalize that result to the case of the arbitrary degeneracy (with 
the Fermi energy ep of any sign and value) it is sufficient to introduce the Fermi distribution 
function in the integrand ((TJ and extend the integration over the interval < k < 00 [S|: 

00 

F( r ,T) = -—[ yN (y)J (y)dy ^ 

r 2 J 1 + exp[(h 2 y 2 /2mr 2 — sp)/kBT] 


The behavior of the function (J2J) is determined not only by the temperature as such but also 
by the temperature dependence of the Fermi energy. The latter circumstance has been ignored 
in j^]. In the framework of the standard two-dimensional band and under the invariable carrier 
concentration, the ratio 77 = ep/k B T (77 < refers to the non-degenerated carriers) is defined 
by the relation 

6 V(T) _ e 7rh 2 n/mkT _ j /g\ 
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that predicts negative 77- values at T = 100 K if n < 10 12 cm -2 . From here, e.g. for the values 
n = 3- 10 12 cm -2 , m = 0.5mo, a = 4A corresponding to the experiments jSlHj, the temperature 
dependence f]{T) results which is shown in the inset of Fig. 1. In fact, it is seen that in the 
real experiment, carriers could not be considered as degenerated (in [Hill] actual temperatures 
amount to 170 K). 

Taking into account (jSJ), the expression could be written in the form 



w(p) = -Je ff <j)(p, t) SiS 2 



where J c s = (ma 2 /4iih 2 )J ( 
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p = r/a is the reduced separation between interacting impurities, r = 2nma 2 kBT/h 2 is the 
reduced temperature 2 , x = N^/N is the relative concentration of magnetic impurities. Under 
the strong degeneracy (n ^> 1), Eq. (JSJ) reproduces the result fTj. while in the case of non- 
degenerate carriers it is similar to the relevant three-dimensional result obtained in In 
the most actual case of the intermediate degeneracy (see below), the function (jHJ) has been 
estimated by numerical calculations. 

In Fig. 1, dependencies 4>(p) are shown calculated with (JSJ) for different values of the carrier 
concentration X7 (the actual range of x7-values is determined by conditions of the concrete 
experiment). Particularly, it is clear that at low enough carrier concentrations (corresponding to 
i] < 0) oscillations of the function <p(p) disappear and its sign corresponds to the ferromagnetic 
interaction. 

For simplicity we use the Ising model corresponding to S = 1/2 and leading, as is known, to 
qualitatively correct results at S > 1, as well. Appropriate generalization does not meet some 
principal difficulties. 

Let the system consisting of randomly arranged and oriented Ising spins be in the state 
characterized by the average reduced magnetization ^ j ^ 1. The total interaction energy 
W = Y2i w i °f a gi ven s P m Si with other spins Sj (i = 2, 3, . . .) is a random value which we 
shall define by the effective local magnetic field H = —W/p (p = gpB\/S(S + 1)) and describe 
by the distribution function F(j;H) depending on the average concentration of effective 
magnetic ions and the reduced system magnetization j = 2£ — 1 where £ is the average fraction 
of spins of magneto-active ions directed up. 

The randomness of the distribution of magneto-active impurities is restricted only by the 
necessity to place them in the fixed locations (sites) of the matrix lattice. For strongly diluted 
systems, that restriction is not significant and the distribution function could be found by the 
Markov's method [T^j, according to which 



F(j;H) 



A(q) exp(-iqH)dq, A(q) = lim 

ZTT I iVmax^OO 



1 N « 

j e iqh ^^K^C)K r {r)dr 
C=±i n 



, (6) 



where h^(r,Q = —C w i r )/p is the field generated at the origin by the spin spaced at the 
random distance r from it. The random parameter ( takes values ±1 (with probabilities £ 



2 For GaMnAs r « 1CT 3 T[K]. 
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and (1 — £), accordingly) and determines the direction of the remote spin, fi^(C) and K r (r) are 
distribution functions for random values of the parameter ( and the distance r, iV max = vrr^^A^ 
is the number of magneto-active impurities in the circle of the radius r max over whose area the 
integration is performed. 

In the spirit of the mean field theory, (^-distribution could be written as 

«c(C) = [(i - 0*(C + 1) + e*(C - 1)]- (7) 

As for the r-distribution, «v(r), one has, in principle, to account for the clustering effect origi- 
nating from the mutual attraction experienced by near Mn-atoms in GaAs-lattice and resulting 
in their non-random (correlated) spatial distribution. Uniform non-correlated r-distribution 
would be 

(2r/(r max — T min ), r > r m j n , , 

r < r min , 



Mr 



where the existence of the minimal separation r m j n between magnetic ions (conditioned by the 
embedding of impurity atoms in the lattice of semiconductor matrix) is accounted for 3 . To take 
into account the correlation of Mn-atoms one has to add the correlation function g{r) in the 
right-hand side of that relation. This function could be, in principle, calculated if the spatial 
dependence of Mn-atoms interaction energy would be known. However, the exact reason for the 
tendency of Mn-atoms to clustering is unclear ^H] ■ At the same time, Monte-Carlo calculations 
[TT] show that impurity correlations have only small effects on the ferromagnetic transition 
temperature of Gai-^Mn^As 3D-system. Thus, though the consideration of the clustering is, 
in principle, straightforward - it is sufficient to introduce in the integrand (JBJ) (and (JOJ), (jTTjl . 
see below) the correlation function g(r) - later on, we will use the simple relation (jSJ) that 
corresponds to g(r) = 1. 

Substituting ©, © in © one finds 



oo 

A(q) — exp[— 2irNfj,C(q)}, C(q) = J {1 — cos[qh(r)] — i ■ j sm[qh(r)]} r dr. 



(9) 



f*min 



Relationships © do not lead to a simple analytical expression for the distribution function 
F x (j; H). So, to determine the latter we have used the low q approximation, based on the fact 
that in the inverse Fourier transform (© the region of high q- values is not important. In that 
approximation 

C(q)=Pq 2 -ijQq, (10) 

where 

oo oo 

P = ~ J h 2 (r)rdr = ( J cS /fi) 2 a 2 0p(p min ), <t>p{r) = \ J ( P 2 (p)P d P^ 

^min Pmin 

(11) 

oo oo 

Q = J h(r)rdr = (J cS /^)a 2 Q (p m i n ), 4>q(t) = J <j>{p)pdp, 

^*min Pmin 



3 The minimal possible distance between magnetic- active Mn-ions substituting for Ga- atoms in the zinc- 
blende AsGa-lattice equals a — a /^/2 w 4 A where clq = 5.7 A is the side of the cubic cell. 
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Ja = 1. 



Substituting (JTUJ), (0) in © we find that in the considered approach the distribution F(j; H) 
is described by the shifted (relative to H = 0) Gauss function 4 



F(j;H) 



2txo 



■ exp 



2a 2 



(12) 



H s = oc N„, a = (2vriV At P) 1 /2 K (13) 

The position of the maximum (H = jHj) of the distribution is determined by the parameter 
Q and depends linearly on the system magnetization j while the distribution width a is defined 
by the parameter P and does not depend on j. The positive sign of Hj means that the 
average direction of the effective magnetic field coincides with the direction of the average 
magnetization, that is the field is, on average, promotes the ferromagnetic ordering of magnetic 
moments. 

Temperature dependencies Hj (r) of the exchange field defined by the relation (JT3Jl are shown 
in Fig. 2. As illustrated in the inset of that Figure, ifj-value peaks at 27 ~ 0.01. It is just this 
condition which results in the maximum Curie temperature (see below). 

Relations (fT3*j) for the shift Hj of the distribution function F(j; H) and its broadening a 
could be rewritten in the form 

Hj = (J eff ///) [27tx 4> q {t)\ , a = (J eir ///) [2nx 0p(t)] 1/2 . (14) 

It follows herefrom 



/J,/a = (2vrx) 1/2 ^(r), V(r) = (15) 



As is shown below, the ferromagnetic ordering is possible under the condition Hj/a > yvr/2 

or 

x 1/2 ^(r) > 1/2. (16) 

Thus, in addition to the material parameters x = N^/N and 7 = n/N^, the function ip{r) 
defining the ratio Hj/a of Gauss function parameters becomes to be crucial in the considered 
problem. As an example, temperature dependencies of that ratio for x — 0.1 and various 7- 
values are shown in Fig. 3. In the inset, the dependencies of the maximum Curie temperature 
r max on ^,_ va j ue following from (J16)) is displayed (for various x- values). 

In the traditional mean field theory, the distribution function is 5-like one for any mag- 
netization j: F(j;H) = S[H — jHj}. It is evident, that the broadening of that distribution 
in a random system prevents ferromagnetic ordering. The magnetization of such a disordered 
system has to be calculated taking into account the scattering of local interaction energies H 
by means of the straight-forward generalization of the equation j = tanh[fiH(j)/kT] referring 
to the regular Ising system: 



tanh 



~kT 



F(j;H)dH (17) 



4 The lattice model ^31 used in [5] corresponds to 5-like distribution function F(j, H) = S(H — jHj) where 
fiHj = Niiufri); Ni, Ti are the number of zth nearest neighbors and their distance. 
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Using the expression ()12|) for the distribution function F(j; H) one gets the equation gen- 
eralizing the standard mean field one: 

-*(?)7-S)- 

— oo 

where 6 = kT//j,Hj = t/xI 2 4>q(t), I = y/irJ ex (h 2 /ma 2 )" 1 is the reduced strength of the inter- 
action 5 . That equation predicts the phase diagram of the system, temperature dependencies of 
its magnetization (in ferromagnetic phase) and susceptibility (in paramagnetic phase), as well 
as the dependence of Curie temperature 8q on the interaction strength, the relative magnetic 
ion concentration x = N^/N and the relative free carrier concentration 7 = n/iV M . 

To clarify under what conditions that equation has the solution corresponding to the ferro- 
magnetic state (j > 0) notice that in the vicinity of Curie temperature where the magnetization 
is small (j —►()), it follows from (JTSj) 

udu = l. (19) 

The integral in peaks at 6 = t = and its maximum value equals (a/Hj) 2 . It follows 
herefrom that the ordered state is only possible under the condition 




which means that the effective RKKY-field proportional to Hj has to overpower not only the 
thermal disordering but also the scattering of local fields proportional to a. The upper boundary 
r max Q f temperature range where the cited condition is satisfied determines the maximum 
attainable temperature of the ferromagnetic ordering at infinite interaction energy (J — > 00). 
It could be derived from the condition 2x 1 / 2 -?/>(t) = 1. The existence of the maximum Curie 
temperature is associated with the lifting of the carrier degeneration at high temperatures and, 
as a consequence, the finiteness of the effective energy of inter-impurity interaction even at 
I — > 00. 

Curie temperature at the finite interaction energy could be determined by solving the equa- 
tion (|19|) . Relevant non-monotone dependencies 70(7) are displayed in Fig. 4. It is clear that 
the relative carrier concentration 7 = n/N^ ambiguously influences Curie temperature in ac- 
cordance with the non-monotone dependence of the exchange field Hj on 7 (cf. Fig. 2). To 
compare, the dashed line in Fig. 4 reproduces the dependence rc(n) presented in and 
obtained in the framework of the standard mean-field theory. 

The optimal carrier concentration occurs to be on the order of 10 12 cm -2 , that is significantly 
lower than the value n ~ 10 14 cm -2 predicted in In addition, Fig. 5 demonstrates there is 
a threshold value of the interaction strength I to drive the system in the ferromagnetic state. 
This is to be contrasted with the result of the standard mean-field theory which predicts no 
such a threshold. 

5 Presently accepted value J p d — 0.15 eV-nm 3 |2] (for GaMnAs) results in / w 1. 
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In conclusion, conditions of establishing the ferromagnetic state and its parameters in quasi- 
two dimensional semiconductor systems with magnetic impurities coupled via RKKY interac- 
tion have been studied in the paper. As distinct from [6 , two new important factors have been 
included in the consideration, allowing for the spatial disarray of interacting magnetic impu- 
rities, and the temperature dependence of the carrier degeneracy. It has been demonstrated 
that both factors complicate transition of the system into the ferromagnetic state: disorder of 
the impurities arrangement reduces the Curie temperature (as compared to the regular system) 
while lifting the degeneracy of carriers makes the Curie temperature finite even in the extreme 
case of the infinitely strong interaction. Besides, the concentration dependence of the transition 
temperature occurs to be non-monotone and there is a threshold interaction strength to drive 
the system in the ferromagnetic state. 

This work has been supported by Grant No. 06-02-116313 of the Russian Foundation of 
Basic Researches. 
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Captions 



Fig. 1. Spatial dependencies 4>(p) of the RKKY interaction energy for magnetic impurities 
in the two-dimensional system with various carrier concentrations degeneracy determining by 
the parameter ary (notice, different scales above and under the Y-axis break). In the insert - 
the temperature dependence rj(T) corresponding to conditions of the experiments [HJE]. 

Fig. 2. Temperature dependencies Hj(t) of effective exchange field for the two-dimensional 
system with various concentrations of magnetic impurities (x = N^/N) and free carriers 
(7 = n/Nfj). In the insert - the maximum attainable if., -value. 

Fig. 3. Temperature dependencies of the ratio Hj/a of Gauss distribution function param- 
eters for two-dimensional system with the concentration of magnetic impurities x = 0.05 at 
various carrier concentrations 7. In the insert - maximum Curie temperature t™ 5 ™ that could 
be attained at given x- and 7-values. 

Fig. 4. Dependencies 7~c(n) of Curie temperature on the carrier concentration for two- 
dimensional system with the concentration of magnetic impurities x = 0.1 for various interac- 
tion strengths I. The dashed line is the result having presented in jHj and corresponding to 

Fig. 5. Dependencies Tc(I) of Curie temperature for two-dimensional system with the con- 
centration of magnetic impurities x = 0.1 at various carrier concentrations determined by the 
parameter 7. 
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